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THEORY OF REVERSIBLE AND NONREVERSIBLE CRACKS IN SOLIDS* 

By Y. I. Frenkel 


THEORY OF GRIFFITH 

The existence of incipient cracks on the surface or inside of solid 
"bodies is, as is known, the reason for the sharp decrease in their actual 
strength in tinilateral tension as compared with the theoretical strength 
value derived from the electron theory of the forces of a molecular chain. 
According to the theory of Griffith (ref. l), in the presence of a crack 
the tensile force applied to the "body, which usually is represented as a 
cylindrical rod, is distributed near the edge of the crack. When this 
’superstress' attains a value corresponding to the theoretical strength 
of the body the crack starts to deepen suddenly as though 'slitting' the 
body; this constitutes the mechanism of its rupture. 

The theory of Griffith is usually presented in a somewhat different 
form, associated with the importance not of the forces but, of the poten- 
tial energy of the body under consideration, more accurately, with the 
change in this energy which is produced by the cracks of given length 
(or depth). In the case of the two-dimensional problem this change, ac- 
cording to Griffith and other authors, may be expressed by a two -termed 
formula 

•pS-j 2 

( 1 ) 

vhere I is the length of the cracky F is the mean value of the tensile 
stress (in the absence of a crack or at a large distance from the latter)^ 
E is Young's modulus, and a is the surface stress. The first term U 
is the decrease in energy due to the formation of cracks of length I (in 
the separating of its edges). The second term represents the correspond- 
ing increase in the surface energy on the assumption that the edges of 
the crack are separated by a distance where the forces of the molecular 
chain may be neglected. The dependence of U on 1 expressed by for- 
mula (l) is graphically represented in figure 1; the mean stress F play- 
ing the part of a parameter. The maximum of the curve U(z) corresponds 
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to the critical value of this parameter =* which the further 

lengthening of the crack hecomes favorable, that is, to the stress failure 
which can be measvired under the usual test conditions. The experimental 
value Z = corresponding to the given value F « F^, characterizes the 

critical length of the crack at which its further lengthening becomes 
favorable. The relation between F]^ and Z]^ is expressed by the well 
known formula of Griffith (derived from the condition Su/9Z = O) 



As increases, Z^ decreases as seen in figure 1; the curves 

F]_ and Fg correspond to two different values of F for the condition 
^2^^1* theory of Griffith the crack is a through hole of ellip- 

tical shape with its major semiaxis equal to Z/2 and its minor semiaxis 
proportional to the mean value of the applied stress F. The magnitude 
of the true stress Fq at the edges of the crack (that is, at the ends of 

the major axis) is inversely proportional to the radius of curvature 6 
of the ellipse at the corresponding vertices and is expressed in terms of 
F]j; by the formula 



For Fq to have a finite value, the radius of curvature 6 must 

likewise be finite. If it is assumed, in conformity with the molec- 
ular concept of the structure of solids, that the minimum value of the 
radius of curvature of the crack should be equal in order of magnitude to 
the distance between the neighboring atoms or molecules (5 = 10“8 cm) and 
if, further, the theoretical strength is identified with the value of 
Young *s modulus E which is in agreement with the electron theory of the 
forces of a molecular chain and if, finally, account is taken that the 
surface stress is equal in order of magnitude to the product of the modu- 
lus E by 5, the value of Fq obtained from formula (3) is practically 

identical with that which is determined by formula (2). 

CRITICISM OF GRIFFITH *S THEORY 

Griffith’s theory is unsatisfactory both in principle and in regard 
to its agreement with experimental facts. Leaving these facts aside for 
the moment, the defects of Griffith’s theory with respect to fundamental 
principle will be considered first in order to construct a theory of 
cracks which is free from defects. 
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It is first necessary to note the following qualification. If for 
a given value F = Fj^ the length of the crack would he found greater 

than the corresponding value (formula (2) or (3)) then, according to 

Griffith’s theory this crack should have increased in a catastrophic 
manner. This conclusion which expresses the essential character of 
Griffith’s theory is generally considered entirely natural and undisput- 
ahle. However, an entirely different result is obtained if the initial 
length of the crack is asstimed less than In this case, according 

to the very meaning of the theory of Griffith, the crack should have spon- 
taneously closed, that is, gradually decreased in length and would have 
finally disappeared entirely. As far as is known to the author, such a 
conclusion from the theory of Griffith has never been drawn by any one 
at any time. Griffith himself assumed, and in this, all his followers 
agreed with him, that for a stress F less than the critical value Fj^, 

which corresponds to the given length of the crack I = the latter 
maintains its length for an unlimited time; in other words, for 7 < 2jj, 

the length I instead of decreasing catastrophically (similarly to its 
catastrophic increase for the case Z > remains constant. This assump- 
tion is a necessary condition for the applicability of Griffith's theory 
to real solid bodies but is at the same time contradictory to it (gener- 
ally without explicit awareness of it). Early in the concept of the 
theory it was established that the initial incipient crack has, even in 
the absence of tensile forces, a finite width considerably greater than 
the interatomic distance 6, representating, as it were, a cut in the body 
not connected with the action of the tensile forces and therefore, incap- 
able of spontaneously closing when the tensile forces are removed. With 
Increase in the latter up to the value Fjj, the Initial crack maintaining 

a constant length can only widen, and only when F exceeds this value 
would it begin to lengthen. 

This contrasting of the initial structure of the crack ('crack with 
hollow’ ) with the structure which it develops on further len^hening is 
however, entirely arbitrary and is Inconsistant with experimental facts. 

In fact, from this supposition the result follows that if on passing 
through the critical value Fj^ and not allowing the crack the possibility 

of starting to lengthen during a certain time interval At, and when the 
applied tensile force is suddenly removed, the crack will spontaneously 
close to its initial length. A test of this kind has not up to this time 
been conducted. There exists, however, every reason to believe that actu- 
ally the initial length of the crack plays no essential part, and that in 
the absence of tensile forces the crack should be entirely closed. 

This defect in the principle of the theory of Griffith may be formal- 
ly reduced to the statement that the curve U(Z) for F = constant has 
only a single extreme point which corresponds to the maximum of the ener- 
gy, that is, to an unstable state of the body in the presence of a given 
tensile force. 
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This qualification deprives Griffith’s theory of any physical mean- 
ing. To obtain a physical meaning of the theory, it is necessary to es- 
tablish a relation between U and Z such that, together with the maxi- 
mum, which corresponds to the unconditional existing limiting value 
Z =s Z", there should be a minimum for a somewhat smaller value Z = Z ' 
which corresponds to the stable state of the body. Such a relation may 
be graphically represented by the ciurve shown in figure 2, which is remi- 
niscent of one of the Isothermals in the theory of Van der Waals (turned 
upside down). 

Curves of precisely this type are, in fact, obtained when the theory 
is constructed on the basis of a true representation of the geometry of 
the cracks and not on the entirely arbitrary and actually incorrect repre- 
sentation by Griffith, which is that elliptical holes with finite curva- 
ture at the ends of the major axis are formed. The present theory related 
to the representation of stable elastic cracks arising from the action of 
tensile forces (external or internal) is in agreement with the investiga- 
tions of P. A. Rebinder on the effect of adsorbing substances or of a 
gaseous medium in which the solid \inder tension is located on its differ- 
ent mechanical properties. This theory corresponds also to the tests of 
Obreimoff on the splitting of sheets of mica along the cleavage planes 
(ref. 2). 


CRACKS FORMED IN AXIAL SPLITTING OF CRYSTALLINE PLATES 

The cracks formed in a solid should not have rounded edges, as is 
implicitly assumed in the theory of Griffith, but sharp edges correspond- 
ing to the gradual asymptotic approach of both surfaces separating the 
crack to the normal distance 6, which in a crystalline body may be de- 
fined as the constant of the crystal lattice in the corresponding direc- 
tion, and which in general is identical in order of magnitude to the 
normal distance between the neighboring atoms or molecules . From a 
purely geometrical point of view these normal distances between neighbor- 
ing layers of particles could be treated as cracks, plane in the case of 
crystalline bodies, and more or less curved in the case of bodies with 
vitreous structure, and the entire solid could be considered as pierced 
in all possible (in particular crystallographic) directions by a system 
of holes of atomic width. However, from a dynamic point of view these 
Intervals between neighboring layers of particles (atoms, molecules) 
become actual cracks only when the distance between them is increased 
relative to the normal distance, that is, becomes equal to S + y; if the 
widening y (the difference between the actual distance and its normal 
value) is considerably greater than 8, even though at a certain limited 
part of the body near the surface or far from it rather than at the corre- 
sponding parts, the edges of the cracks may be considered completely dis- 
connected. Such a slit-shaped crack is represented in figure 3(a) for the 
case of a crack on the surface of the body and in figure 3(b) for the case 
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of an internal crack. The abscissa here denotes the distance measured 
along the crack and the ordinate the corresponding widening y or, more 
acciirately, half of it (thus the value y =« 0 characterizes the normal 
distance hetween the nei^horing layers of particles ).l 

Although under the observed conditions, the length of the crack 
strictly spe a king is infinite, practically it may be assumed that the 
crack, in the actual sense of the word, ends where the distance between 
its edges becomes equal, say, to twice the value of the normal distance 
6, that is, where y = 5. 

A less arbitrary and physically more justified definition of the 
effective length of the crack may be obtained on the basis of the depend- 
ence of the mutual interaction between the oppositely lying parts of the 
walls of the crack, and on the distance between them y (or, more accu- 
rately, 6 + y). This force, referred to unit area, is denoted by f(y). 
For sufficiently small (in comparison with 6) values of y, the force 
should, evidently, be proportional to y, reducing to an attraction in 
the case y>0; and to a repulsion in the case y<0. This dependence 
may be expressed by the formula: 


f 



( 4 ) 


where f >0 corresponds to an attraction, and f<0 to a repulsion. 
Since for y > > 6, f is equal to 0, it follows that for a certain 
value y = y* the force f attains its maxinnim value f*, which may be 
identified with the theoretical strength of the material in rupture. 

Since y* is near in value to 5, f* should be near in value to E, as 
is clear from formula (4) which may be approximately extended to this 
case but becomes entirely unsuitable for y>y* 

The dependence of f on y is graphically shown in figure 4. It 
may be analytically represented by the very simple formula 

f = Aye“^y (5) 

where A = e/S and the coefficient a is equal to the reciprocal of the 
distance y*. For simplicity in this report a will be identified with 
1 / 6 . 


^he difference "between our proposed description of the crack and 
the description that follows from the theory of Griffith does not require 
explanation. We may remark only that the ‘rounding’ of the edges of the 
Griffith cracks with a radius of curvature 5 would correspond, from 
our point of view, to the cracks for which the widening y, instead of 
asymptotically approaching zero, becomes zero in a more or less sharp 
manner. The assuirption made by equating the minimum value of the radius 
of curvature to 6 is entirely without foundation. 
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Formula (5) may "be replaced "by a two-term expression of the type 

f = 5 — (v - |a > 0) (5a) 

(6 + y)^ (B + y) 

in which the first term characterizes the forces of attraction and the 
second (generally omitted in considering the surface energy), the forces 
of repulsion. 

Evidently^ integrating f with respect to y between the limits 
0 and results in the work of the total division of the body into two 
parts, that is, twice the value of the surface energy (or of the surface 
stress in the case of absolute zero temperature). Therefore, 

If, as the upper limit in this integral an infinite value of y is 
taken as a measure of the surface energy for the cracks (slits) not having 
entirely separated edges, this expression applies 

ScrCy^^) = fdy 


As in the case of formula (5) 


2o^(yn) = -o 1 « (ynOL + l)e 


For this expression reduces to 20 q = A/a^; hence, in the 

relation A = E/5 (that is. A/a = E), there is obtained the formula 

cTo = i E6 (6a) 

which practically agrees with the relation between the surface stress and 
the rupture strength that was identified earlier with E. 


More accurately, this strength is expressed by the formula 

f* = - e = Ee = 2.72E 
a 


The accuracy is, however, illusory because formula (5) possesses a roughly 
approximate character the same as the identification of y*(= l/a) with 6. 

When the slit-shaped crack represented in figure 3(a) or (b) is 
taken into consideration, it is assumed that the surface energy of its 
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principal part (corresponding to y>5) is equal in unit width to SuqZ, 

where I is the effective length; while the surface energy of the deeper 
part of the crack may he represented with a sufficient degree of accuracy 
hy the integral of the expression 



taken along x from x « Z to x ** •: 

'’-ill' 

That the siirface energy of a narrow crack should be less than that 
of a wide crack was noted by Rebinder in his work on the strength of 
solids under the action of substances adsorbed on their surfaces. 
Rebinder, however, restricted himself to only a qualitative determination 
of the function a(y) . 


DEPENDENCE OF LENGTH OF CRACK ON THE MAGNITUDE OF 
TENSILE FORCES FOR A SPLIT PLATE 

Imagine a horizontal plate of infinite length split under the effect 
of two oppositely directed forces iF applied perpendicular to its plane 
and to its edge from above and from below, respectively (fig. 5). 

The width of the plate (in the direction perpendicular to the plane) is 
taken to be unity. The thickness of the plate is denoted by 2b which 
is assumed to be small in comparison with the effective length I of 
the forming crack (lengthened in the form of a gradually narrowing slit). 

To determine I as a function of F recjuires that each of the two 
similar halves into which the plate splits be considered as a beam bending 
under the action of a force concentrated on the free end and the cohesive 
forces on the other half of the plate, which are distributed along the 
surface directed toward this half, and attaining a maximum value at the 
required distance 2. In the work of Obreimoff (ref. 2), account was taken 
only of the wide part of the crack or its *mouth*, while its slit-shaped 
part, extending from x =* 2 to x = •, which is the usual concept of a 
crack of finite length, was not taken into account at all. These surface 
forces take the place of the forces which in the usual theory of the bend- 
ing of a beam are reduced to the reaction of the wall to which one of its 
ends is clamped. 

Under these conditions the shape of the plate, or more accurately, 
of one of its halves (the lower for example) is determined by the differ- 
ential equation 
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( 8 ) 


where y denotes the displacement of the corresponding point of the 
half -plate under consideration^ "because the lower half -plate is taken to 
he positive downwards^ E is Young's modulus^ I = h^/3 is the moment of 
inertia of a plane section^ and M the moment of the forces applied to 
the half -plate with respect to the point x = 0. This moment is expressed 
hy the following formula 


M = Fx 



f 4 ^ 


(8a) 


where f is the surface force referred to unit length. Differentiating 
this expression with respect to x gives 


m 

dx 


= F 


- fx 


(8b) 


Equation (s) is thus^ equivalent to the following third-order 
equation 


d^y 


= ^ (F - fx) 


( 9 ) 


where f is a function of y determined hy formula (5) or (5a). 


For the solution to the problem of the dependence of the effective 
length of the crack on the external force F^ an exact solution of equa- 
tion ( 9 ) is not necessary. The problem may he solved hy an approximate 
method taking into account that for large values of x the derivative 
d^y/dx^ should approach zero. Since it is then possible to put 


E F 

f=Ay = ^y=— ^ the dependence of y on x in the slit -shaped part 

O X 

should be expressible by the asyiiiptotic formula 


F S 
^ Ex 


( 10 ) 


obtained from equation (9) under the condition 
from equation (lO) it follows that d^/dx^ = - 


the following corrected expression 

F5 



d^y/ dx^ 
6F5 2^ 
E ^4^ 


= 0. In fact^ 
which leads to 


(10a) 


For the purpose of this report the correction term will he neglected. 
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It is understood "by the effective length I of the crack that the 
value of X, for which the value of y is equal to 5, corresponds (ap- 
proximately) to the maximum value of the force. Thus the following linear 
-relation between I and F can be obtained 


I 



(11) 


ENERGY OF CRACK IN SPLITTING OF A PLATE AND STRENGTH OF THE LATTER 


The preceding results may also be obtained from the condition of the 
maximum total energy of the system under consideration (half -plate ) . This 
energy U breaks up into three parts; (a) the elastic energy of the 

external part treated as a beam of length I with load F on the free 
end; (b) surface energy of the open part of the crack Ug = and (c) 

the energy of the incipient slit-shaped part corresponding to x >2. For 
part (c), formula (7) may be used with formula (lO), which gives 


U 


3 



( 12 ) 


Part (a) is obtained by the expression 

Ui = - I Fyo 


where y^ is the maximum deflection, that is, the value of y for 
X =s 0. This may be expressed by the approximate formula 


^0 


Z^F 

3EI 


which is obtained by the integration of equation (s) substituting in it 
M ss Fx, that is, not taking into account explicity the internal forces 
but schematically using the conditions y ** 0 and dy/dx » 0 for x « Z. 

Thus, the total energy of the plate, or more accurately, the half- 
plate of the splitting force F over the effective depth Z may be ap- 
proximately represented by the following formula 


U 


F^Z^ 


20^1 + 


F^5 

EZ 


(13) 


which of coTirse is applicable to the case of very small values of e 
(These values will not, in this report be taken into account.) 
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The dependence of U on Z expressed by formula (l3) is graphi- 
cally shown by the curve of the same type as in figure 3 with the differ- 
ence that when Z approaches 0, U approaches • (in reality it ap- 
proaches a finite value). In figure 6 is shown a family of curves U(Z) 
corresponding to different values of the splitting force F (isodynes) 
in the order of increase of the latter. For a value that is not too large 
the curve has two extremities; a minimum for Z = Z' and a maximum for 
Z = Z" > Z ' . By increasing F these points approach each other and 
finally for a certain value F * F^. (the critical point analogous to the 

critical temperature in the theory of Van der Waals) they coalesce. 

For F<<Fg, the miniimim of the function U practically coincides 
with the minimum Of the sum 

F^& 

Uj + U3 = 20qI + 

and corresponds to the following value 



Since, according to equation (6a) cTq =* E8/2 this expression may be 
written in the form 

Z* = I (14a) 

which agrees with the previously derived expression (ll) for the effective 
length of the crack due to the action of the force F. From this it is 
clear that the length of cracky which was considered earlier^ represents 
a stable value, corresponding to the miniinum of the potential energy U 
and characterized hy the stable equilibrium of the split plate. 

The maximum value of U corresponds to the unstable equilibrium of 
the plate at which the length of the crack, having attained a certain 
value 7" > 7 ’ under the action of a very small distiorbance, either re- 
turns to the value I =7*, or suddenly increases. For the condition 
F 5= F^, this length practically reduces to that which corresponds to the 

3 

maximum of the function U]_ + U 2 + 2aQ7 and is expressed by the 

formula 

= 1 = 


Z" 


(15) 
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Corapaxing equations (l4a) aud (15) we may write the relation between I 
and Z* in the form 


Z'Z” 



(15a) 


Evidently the value corresponding to the agreement of Z • and 

Z” simply represents the strength limit of the plate relative to the 
action considered, that is, the minimum force introduced by the sudden 
formation of the lengthening crack. ^ In order of magnitude, this force 
may be determined by identifying expression (l4a) with (15) . Still more 
simply it may be computed with the aid of formula 

^c=VF^ (16) 


and 


F 


c 


EZc 


(16a) 


When 6 
meter and F^ 


10"6 centimeter and b = 1 centimeter, Z» *« 10“^ centi- 
+10 , +12 , 2 > ° 

10 dynes (for E = 10 dynes/cm ) are obtained. 
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In contrast to the theory of Griffith this force is not connected 
with the existence of an already existing crack of a certain finite 
length. 
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